© 

NEL  Memonuidam  Report  5194 

Additional  Capabilities  of  the  Unified  Waves  Model 

J.  M.  Witting 

Laboratory  for  Computational  Physics 


September  29,  1983 


NAVAL  RESEARCH  LABORATORY 
WashiogtOD,  D.C. 

Approved  for  puUk  reieese;  distribution  unlimited. 

94  2  24  180 


AD-A276  114 

llililllllil 


Best 

Available 

Copy 


SCCUMlITV  CLikJSiriCATlON  OF  This  PkUC  CWliwi  0«a  Inftmai 


REPORT  DOCUMENTATION  PAGE 


MC^OMT  NUMACA 


;a  COVT  ACCESSION  no 


NRL  Memorandum  Report  5194 


4  TITI.C  (in 

ADDITIONAL  CAPABILITIES  OF  THE  UNIFIED 
WAVES  MODEL 


7.  AU  TMOm(t/ 


J.  M.  WitUng 


KEAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1  RECIPIENT'S  catalog  NUMRER 


S.  TYPE  OF  REPORT  •  PCRlOO  COVEREO 

Interim  report  on  a  continuing 
NRL  problem. 


f  performing  ORG.  REPORT  NUMRER 


*  contract  or  grant  nuMRERT*) 


t.  performing  organization  name  ANO  AOORESS 


Naval  Research  Laboratory 
Washington,  DC  20375 


10.  PROGRAM  element,  project.  TASK 
AREA  •  RORK  UNIT  NUMRERS 

61153N;  RR  032-04-03; 
44-1868-03 


II.  CONTROLLING  OFFICE  NAME  AnO  AOORESS 


OfRce  of  Naval  Research 
Arlington,  VA  22217 


monitoring  agency  name  a  AOORESSrK  RlffaraAl  I 


I  C«ni»o/llii4  OtHc*) 


12.  NC^ONT  DATS 

September  29, 1983 


n  numbcn  of  paces 


IS.  SeCUNlTV  CLASS.  f»i  (Af# 

UNCLASSIFIED 


If-  OlSTNieutlON  STATCMCNT  (of  (Mi  JI«bap(> 

Approved  for  public  release;  distribution  unlimited. 


17.  OlSTRieuTiON  statement  (Of  ttto  «di(r«<l  ontoroii  in  Block  20.  tt  dttloront  from  Rop^rt) 


If.  KEY  WOROS  fC«n(litw#  on  rooofio  oiOc  II  nooooiovy  IBontliy  ky  AlocA  wuiAor^ 

Waterwaves  Steady  nonuniform  flow 

Unified  waves  model 

Nonlinear 

Variable  depth 

Wavemaker 


20.  AaSTRACT  fConllnwo  on  rovoroo  olBo  II  noooooorr  on#  IBonilty  kf  AlocA  nuoiAora 

The  initial  report  on  the  “uniOed  waves  model”  describes  a  procedure  to  compute  water  wave 
propagation  over  long  distances  accurately,  at  tow  cost,  and  for  a  wide  variety  of  physical  situations 
(Witting,  1982).  Iliis  report  outlines  recent  developments  that  extend  its  capabilities.  These  involve: 
(1)  inclusion  of  a  variable  depth,  (2)  addition  of  a  wavemaker,  and  (3)  the  ability  to  send  steady  non- 
uniform  flows  through  the  computational  domain.  Examples  of  calculations  employing  the  new 
capabilities  are  given. 


00  /A 


EDITION  OF  I  NOV  ES  IS  OaSOLETe 
S/N  OIOZ-014-  4601 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (Whit  OM.  RnlwMI 


I 


CONTENTS 


1.  INTRODUCTION  .  1 

2.  DISPLAYS  OF  THE  NEW  CAPABIUTIES  .  1 

3.  DISCUSSION  OF  THE  NEW  CAPABIUTIES  .  5 

REFERENCES  .  9 

ACKNOWLEDGMENTS .  9 


1  Accesion  for  j  j 

NTIS 

CRA&I 

^ - 

OTIC 

TAB 

□ 

j  Unannoiir'.ced 

□ 

1  JiJStifiCJt  ii* 

'ev 

Distribulion  / 
_ 

A 

v.Mi.'h'htv  Tofles  1 

i  Awdll 

jiid/or  1 

Oist 

'  Coecial 

# 

iii 


ADDITIONAL  CAPABCLITIES  OP  THE  UNIFIED  WAVES  MODEL 


1.  Introduction 

This  brief  report  outlines  developments  that  have  been  added  to  the 
unified  waves  model  to  extend  its  capability  beyond  that  described  by 
Witting  (1982).  These  include:  1)  inclusion  of  a  variable  depth, 

2)  addition  of  a  wavemaker,  and  3)  inflow-outflow  boundary  conditions. 
Section  2  displays  the  new  capabilities,  and  Section  3  provides  details  of 
the  developments. 

2.  Displays  of  the  New  Capabilities 

Figure  1  shows  results  for  a  run  in  a  channel  of  variable  depth,  with  a 
wave  train  generated  by  a  wavemaker.  Between  x  ■  80  and  x  *  160  the  channel 
shoals  linearly  from  one  depth  to  a  depth  ten  times  smaller.  Distances  are 
non-dimensionalized  with  the  deep  undisturbed  depth  being  unity  (Fig.  la). 
The  water  is  initially  undisturbed  (See  lowermost  trace  of  Fig.  lb).  A 
wavemaker  is  then  turned  on  to  produce  a  moderate-amplitude  wave  traveling 
away  from  the  rigid  left-hand  wall.  Figure  lb  is  a  coarse  (every  fourth 
grid  point)  snapshot  of  the  wave  elevation  at  the  initial  and  32  equally- 
spaced  time  levels.  The  spacing  is  such  that  a  disturbance  that  travels  at 
a  speedy^h^  (h  is  the  depth  at  the  left)  would  show  up  along  a  ±  45° 
angle  in  the  figure.  When  the  wave  shoals  it  gets  higher  and  thinner. 

Before  the  end  of  the  run  the  numerical  results  show  a  wave  with  "higher 
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than  breaking"  slopes.  The  code  crashed  four  time  steps  after  the  uppermost 
profile  .* 

The  remainder  of  Fig.  1  gives  details.  Figures  (Ic-lj)  show  that  the 
wavemaker  is  working  satisfactorily.  These  are  blown-up  views  of  the  5th, 
9th...,  33rd  profiles  of  Fig  lb,  for  the  leftmost  600  grid  points  (there  are 
1600  grid  points  in  the  conputational  box;  Ax  »  At  =  1/8  for  this 
calculation).  The  wavemaker  produces  smooth  waves.  Details  of  the 
wavemaker  will  be  given  in  Section  38. 

Figures  (Ik-lr)  center  attention  on  the  next  600  grid  points.  The 
shoaling  starts  at  grid  point  640  and  continues  to  grid  point  1280.  The 
29th  profile  (Fig.  Iq)  is  used  for  discussion.  It  is  a  continuation  of  Fig. 
li.  The  leftmost  wave  in  Fig.  Iq  is  a  continuation  of  the  steadily 
marching  waves  in  Fig.  li.  The  effects  of  shoaling  are  apparent  in  the  rest 
of  Fig.  Iq.  The  waves  get  higher,  and  travel  slower.  They  also  become 
sharp-crested. 


♦The  analysis  and  numerical  methods  conserve  exactly  those  linear  quantities 
that  should  be  conserved  from  the  exact  prognostic  equations  resulting  from 
applying  the  continuity  and  Euler  equations  to  the  situation  (Witting  and 
McDonald,  1982  derive  the  form  used).  These  conserved  quantities  are  the 
integrated  mass  and  an  integrated  surface  velocity.  The  model  involves 
approximations,  however,  and  physical  energy  is  not  conserved  exactly, 
though  the  code  is  reversible,  so  that  no  "numerical"  dissipation  is 
introduced.  The  experience  gained  from  running  the  model  in  a  variety  of 
problems  is  that  energy  is  conserved  rather  precisely  for  physical 
conditions  where  waves  are  not  breaking,  but  not  so  precisely  for  physical 
conditions  where  the  waves  would  be  breaking.  Indeed,  the  model  usually  .  ' 
becomes  rapidly  unstable  at  wave  amplitudes  near  (usually  above)  the 
breaking  point.  The  inclusion  of  local  dissipation  to  simulate  breaking 
waves  would  be  desireable  both  to  improve  the  model  and  to  stabilize  it. 
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Figure  Ir  is  taken  36  time  steps  later.  The  wave  between  1050  and  1100 
is  very  sharp-crested  indeed.  Four  time  steps  later  the  instability  that  is 
seen  starting  here  grew  to  floating-point  overflow  catastrophe. 

Figures  (Is-lz)  concentrate  on  the  uppermost  8  profiles  of  Figure  lb. 
Note  changes  in  both  vertical  and  horizontal  scale  from  the  previous 
figures.  The  last  two  (Figs,  ly-lz)  show  the  sharp-crestedness. 

Figures  (laa-lee)  center  attention  near  grid  point  1075,  where  the  wave 

^  ' 

breaking  eventually  occurs,  for  the  uppermost  five  profiles  of  Fig.  lb.  The 
scales  are  set  so  that  the  profiles  are  not  vertically  exaggerated,  and  the 
location  of  the  shoaling  bottom  is  shown.  As  we  proceed  from  Fig.  laa 
through  Fig.  Idd  note  the  solitary-wave-like  appearance,  with  increasing 
amplitude.  The  aoplitude  and  maximum  slopes  of  Fig.  Idd  exceed  those  for 
breaking  progressive  waves.  Fig,  lee  shows  what  is  about  to  happen  to  the 
numerics  when  the  model  wave  height  is  about  twice  that  of  the  highest  wave 
in  water.  Details  of  how  the  variable  depth  has  been  added  to  the  unified 
waves  model  are  described  in  Section  3A. 

Figure  2a  is  a  display  of  the  inflow-outflow  capability.  A  steady  flow 
is  directed  from  right  to  left  in  a  channel  of  variable  breadth  that  is 
narrower  at  the  right  than  at  the  left.  Consequently,  the  flow  is 
decelerating  from  right  to  left.  A  wave  packet  initially  centered  at  x*50, 
is  designed  to  move  to  the  right  (but  a  little  disturbance  also  moves  to  the 
left).  Figure  2a  shows  the  resulting  fate  of  the  packet  (its  caption 
sunoiarizes  what  happens).  As  can  be  seen  from  the  figure,  the  waves  are 
blocked  by  the  adverse  currents.  Linear  theory  predicts  the  blocking 
location  to  be  at  x*119. 
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The  calculations  for  Fig.  2  differ  from  those  of  Fig.  I  and  all  those  in 
Witting  (1982)  in  that  a  "trapezoidal  correction"is  added  at  each  time  step. 
This  procedure  is  slightly  dissipative,  suppressing  an  alternating  grid 
point  instability  that  is  bothersome  when  steady  currents  run  through  the 
computational  channel.  Consequently,  it  is  not  possible  to  track  the  wave 
packet's  energy  with  the  data  shown  in  Fig.  2.  A  successful  effort  has  been 
made  in  subsequent  calculations  that  suppress  instabilities  toward  the  right 
of  x*150  but  are  dissipationless  for  x<150,  over  which  all  the  action  takes 
place.  These  run  to  longer  times.  Results  agree  with  those  of  Figure  2. 
They  show  that  the  packet  reflects  no  visible  waves,  and  remains  nearly 
stationary  with  its  leading  edge  near  x»120  (phases  enter  at  the  left  and 
disappear  at  the  right  of  this  dispersive  wave  train).  Ao^litudes  and 
slopes  at  the  front  of  the  train  eventually  exceed  those  of  breaking  waves, 
and  so  the  physical  situation  would  be  that  of  deep-water  whitecapping.  It 
may  be  concluded  that  the  fate  of  a  wave  packet  encountering  sufficiently 
strong  adverse  currents  is  to  break. 

Figures  (2b-2i)  are  centered  at  the  point  at  which  the  steady  current 
gradients  are  near  maximum.  On  display  are  blowups  of  the  5th,  9th, ...  ,33rd 
(uppermost)  curve  of  Fig.  2a.  Note  the  decrease  in  wavelength,  the  increase 
in  amplitude,  and  eventual  sharp-crestedness  of  the  waves  in  the  center  of 
the  train.  Some  of  the  increase  in  amplitude  is  caused  by  the  waves  passing 
through  a  gradually  narrower  channel.  The  majority  is  due  to  wave-current 
interaction.  Details  of  how  the  inflow-outflow  capability  has  been  added  to 
the  model  are  given  in  Section  3C. 
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3.  Pi»cu8»ion  of  the  New  Capabilities 
3A.  Variable  Depth 


The  analysis  presented  in  Witting  (1982)  is  valid  for  waves  in  channels 
of  gradually  varying  breadth  and  depth.  In  the  initial  computer 
implementation,  however,  the  depth  was  set  to  unity,  and  the  "h"  that 
appeared  in  the  analysis  was  replaced  by  l+n(x, t)  To  get  the  depth  into  the 
problem,  it  was  necessary  only  to  replace  l  +  il(x,t)  with  h^Cx)  +  n(x,t), 
where  h^Cx)  is  the  undisturbed  depth  normalized  to  a  depth  somewhere  in 
the  channel,  x^,  say,  where  h^(x^)*l.  The  undisturbed  depths  h^Cx) 
are  supplied  to  the  computer  externally  (by  assigning  a  file  and  reading 
it).  No  dif f iculities  were  experienced  in  the  course  of  adding  the  variable 
depth  capability. 

The  hQ(x)  used  in  Fig.  I  is  a  set  of  three  straight  lines  with 
intersections  smoothed  over  a  distance  of  approprimately  0.5  hQ(xQ).  No 
tests  were  performed  to  see  whether  unsmoothed  variable  depth  channels  would 
pose  numerical  problems. 

3B .  Wavemaker 

A  wavemaker  was  needed  to  be  able  to  continuously  deposit  energy  with 
known  frequency  and  wavelength  content  into  the  channel,  preferably  near  one 
end.  Pumping  water  into  and  out  of  the  channel  from  an  end  produced  waves 
with  ragged  waveforms.  An  alternate  procedure,  punping  water  without 
horizontal  momentum  into  and  out  of  the  channel  at  several  grid  points  near 
an  end  produces  waveforms  of  the  quality  of  those  shown  in  Figs  (Ic-lj). 

The  procedure  for  designing  the  wavemaker  employs  linear  water  wave 
theory: 

(1) 
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where  n  is  the  elevation  above  still  water  level,  u  is  the  vertically 


(2) 


averaged  horizontal  velocity,  is  the  horizontal  component  of  surface 
velocity,  Q(x,t)  is  a  mass  source  (the  zero  on  the  right  side  of  (2)  inplies' 
that  no  horizontal  momentum  is  injected  into  the  fluid).  For  linear 
sinusoidal  waves  one  can  show  that 


u  =  u  tanh  (kh  )/kh  . 
s  o  o 


(3) 


We  want  to  find  Q(x, t)  such  that  a  nice  wave  comes  out.  We  do  this  by 
writing  the  elevation  as  a  function  that  behaves  like  a  progressive  wave  as 
x+«*  and  as  a  standing  wave  with  antinode  at  x*0  at  the  wall  where  the 
velocity  vanishes.  We  take  a  steady-state  solution.  For  Fig.  1  the  form 
chosen  is: 


(cos  ut  cos  kx  +  sin  mt  sin  kx  tanh  one} 


(4) 


The  next  step  is  to  evaluate  u^  using  (4)  and  (2).  Then  evaluate  u 
from  Ug  using  (3) . 

Finally,  after  algebra,  (1)  gives: 

Q  =  2<A)n  cos  ojt  sech^Qcx  cos  kx  -  r~r  sin  kx  tanh  ok}  (5) 

^  o  k  k 

In  the  numerical  model  source  terms  are  inserted  to  permit  mass 
injection  at  each  time  step  of  the  form  of  (5),  with  k,  and  a  beig 
inputted  parameters,  and  u  being  evaluated  from  the  linear  dispersion 
relation. 

The  data  used  to  make  the  wave  shown  in  Fig.  1  are: 

n 

kh  *  1.0;  =  0.1;  oh  =  0.1  (6) 

o  h  o 

o 
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After  a  abort  time  the  wave  train  haa  acquired  aome  of  the 


characteristica  of  (4)  in  its  main  body  (note  that  (4)  is  a  steady-state 
solution,  and  our  waveraaker  has  for  time  <0,  and  Q  from  Eq.  (5)  for 
t>0).  In  particular,  the  wave  train  has  the  expected  wavelength  and 
frequency.  Its  amplitude  away  from  the  wall  exceeds  that  expected  by  tens 
of  percent,  however.  This  may  be  due  to  applying  a  linear,  steady  state 
analysis  for  waves  that  are  nonlinear  and  initially  transient.  In  any  case, 
the  wavemaker  of  (5)  produces  acceptably  clean  progressive  waves  of  known 
frequency-wavelength  content,  which  is  good  enough. 

3C.  Inflow-Outflow  Conditions. 

The  goal  of  this  development  is  to  allow  a  steady  current  to  flow 
through  the  conputational  channel;  this  differs  from  a  frequently-required 
goal  in  computations,  to  let  disturbances  radiate  out.  The  only  coding 
change  necessary  in  the  running  programs  is  to  ensure  that  end  wall  boundary 
conditions  are  set  to  give  a  correct  mass  transport,  surface  velocity,  and 
elevation  that  are  appropriate  for  a  steady,  undisturbed  flow.  Consistency 
requires  that  if 

q  =  u  +  n'v  =  u  =  u  ,  (7) 

^s  s  s  s  o 

<|»  =  u  b  dy  =  u  bh  =  u  b  [h  +n  (t)  )  (8) 

so  ^  o  o  o  o 

at  the  wall  (or  vice  versa). 


Providing  suitable  input  files  is  a  little  more  subtle.  Assuming  that 
volume  fluxes  into  and  out  of  the  channel  are  the  same,  we  can  have  a  steady 
flow  in  a  given  geometry,  but  only  if  the  initial  surface  elevations  and 
velocities  are  set  at  precise  values;  otherwise  long  waves  propagate.  These 
long  waves  would  bounce  around  the  channel  at  speeds  of  the  order  of  unity, 
and  upset  both  numerical  stability  and  interpretation.  After  some  trial  and 
error,  a  satisfactory  way  of  handling  the  input  files  was  found.  A  channel 
geosietry  can  be  selected  and  long  wave  theory  used  to  give  n(x,  t»0)  that 
corresponds  to  a  steady  flow  in  this  theory.  From  n(x,t“0)  and  the  known 
velocity  at  one  wall  one  can  compute  exactly  q_  (x,t«0)  for  a  steady  flow 
from  Bernoulli's  Law.  From  n  and  one  then  can  recompute  the  channel 
geometry.  The  new  geometry  should  differ  only  slightly  from  the  original, 
but  should  keep  the  long  wave  steady. 

The  above  procedure  was  used  in  part,  for  a  channel  of  variable  breadth 
(I  have  not  yet  tried  it  for  a  channel  of  variable  depth).  A  smooth  surface 
elevation  was  assumed  for  the  steady  flow  of  the  form: 

n(N)  -  -0.25  (An)  tanh^  (8(H-800)  +ll  (9) 

where  N  is  the  grid  point  number.  An  is  the  magnitude  of  the  depth  change 
from  end  to  end  (*0.075  for  Fig.  2),  and  6  expresses  how  gradually  the 
change  occurs  ("0.005  for  Fig.  2).  The  velocity  at  the  left  was  set  to 
-0.100.  With  the  profile  of  (9)  it  went  to  -0.400  at  the  right.  The 
channel  breadth  that  permits  this  is  smoothly  varying,  between  b“1.0  taken 
at  the  left  to  b"0.270  at  the  right.  This  channel  breadth  configuration  and 
these  steady  flow  conditions  are  shown  in  Fig.  3. 


Onto  Che  steady  flow  can  be  added  either  a  waveraaker  or  a  superposed 
initial  disturbance  well  away  from  walls.  The  latter  procedure  was  used  for 
Fig.  2. 

One  may  inquire  what  physical  conditions  match  the  endwall  boundary 
conditions  of  the  model.  It  appears  that  they  are  rigid  walls  for  high- 
frequency  waves,  yet  permit  steady  flow  through  them.  A  poraus,  dense  wall 
with  high  constant  pressure  head  driving  Che  steady  flow  seems  Co  be  a 
possible  physical  condition  chat  corresponds  this  numerical  implementation 
of  Che  unified  waves  model. 
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Fig.  la  and  b  —  Wavaa  ganaiatad  bjr  a  wavemakar  propagating  through  a  channel  of  decreasing  depth ; 

(a)  the  geometry  and  (b)  eucceaeiee  profilea  of  eleeation,  etarting  at  the  bottom.  A  wavemaker  is  turned 
on  at  t "  0  and  makes  meas  at  a  fixed  frequency  each  that  the  initial  wavelength  is  2irho,  hg  being  the 
undisturbed  depth  at  the  left.  Hie  spacing  between  sueceasiTe  profiles  is  scaled  so  that  a  feature  that 
would  travel  at  a  speed  Vgho  would  appear  along  a  ±46*  line  on  the  flgure.  The  calculations  are  run 
widi  Ax  ■  At  *  1/8  for  1792  time  stops  (to  t  ■  224)  in  a  computational  region  of  1600  grid  points 
(to  X-  200). 
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Fig.  1(c)  DetaiU  of  the  scene  near  the  left  well  at  t  «  28 
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Fig.  1(d)  ^  Details  of  the  scene  near  the  left  wall  at  t «  56 
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1(S)  ~  Detaila  of  the  acene  near  the  left  wall  at  t  -  140 
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Fig.  1(h)  —  Detaila  of  the  acene  near  the  left  wall  at  t  •  168 
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Fig.  l(i)  —  Detaila  of  the  scene  near  the  left  wall  at  t  -  196 
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Fig.  l(j)  —  Details  of  the  scrae  near  the  left  wall  at  t  >  224 
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ng.  l(m)  —  Detaila  of  the  •ecne  over  the  ihoaling  portion  of  the  channel,  between  z  ■  76 
(grid  point  600)  and  x  ■  160  (grid  point  1200).  Tlie  ihoaling  itarta  at  z  -  80  (grid  point 
640)  and  oontinuei  to  z  *  160  (grid  point  1280,  outride  theie  figurei)  at  t  -  84 
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Fig.  l(n)  —  Detail!  of  the  scene  over  the  shoaling  portion  of  the  channel,  between  z  »  76 
(grid  point  600)  and  z  -  160  (grid  point  1200).  Tlie  shoaling  starts  at  z  -  80  (grid  point 
640)  and  continues  to  z  <■  160  (grid  point  128>0,  outside  these  flgures)at  t  >-112 


l(<l)  ~  Detaila  of  tha  scena  over  the  ahoaiing  portion  of  the  channel,  between  a  •  75 
(grid  point  600)  and  x  >  160  (grid  point  1200).  The  ahoaling  atarta  at  a  >  80  (grid  point 
640)  and  continuea  to  a  ■  160  (grid  point  1280,  outaide  theae  figurea)  at  t  >  196 
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Fig.  l(r)  —  Oetaila  of  the  acene  over  the  ahoaling  portion  of  the  channel,  between  a  >  76 
(grid  point  600)  and  a  -  160  (grid  point  1200).  The  ahoaling  atarta  at  a  -  80  (grid  point 
640)  and  eoathraea  to  a  ••  160  (grid  point  1280,  outaide  theae  figurea)  at  t  >  224 
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F1(.  1(f)  —  Further  detailf  of  the  aceae  between  x  ■  118.76  (|rid  point  960) 
and  X  <■  166.26  (grid  point  1260)  at  t  -  176 
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Fig.  l(t)  —  Further  detaib  of  the  aeene  between  x  118.76  (grid  point  960) 
and  X  -  166.28  (grid  point  1260)  at  t  -  182 
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Fig.  l(w)  —  Further  detail*  of  the  icene  between  *  •  118.76  (grid  point  960) 
and  X  -  156.25  (grid  point  1250)  at  t  =  203 


Fig  l(x)  —  Further  detail*  of  the  *cen*  between  x  -  118.76  (grid  point  950) 
and  X  -  1 66.26  (grid  point  1 250)  at  t  -  210 


SHIO  rOINt 

Pic*  l(y) Purthar  daUilt  of  tha  acana  batwaan  x  ■  118.76  (grid  point  960) 
and  X  -  166.26  (grid  point  1260)  at  t  -  217 
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Fig.  l(s)  —  Furthar  dataila  of  tha  acana  batwaan  x  -  118.76  (grid  point  960) 
and  X  -  166.26  (grid  point  1260)  at  t  ^  224 
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Fic.  l(aa)  —  The  scene  between  x  ■  132.376  (grid  point  1063)  and  x  <•  135.76  (grid  point  1086). 
The  horizontal  and  vertical  scales  are  the  same;  at  t  »  196. 


SRIO  fOINT 

Fig.  l(dd)  —  Tile  ecene  between  x  •  132.375  (grid  point  1063)  and  z  -  135.75  (grid  point  1086). 

Hie  horizontal  and  vertical  acales  are  the  same:  at  t  •  217.  The  wave’s  maximum  slope  and  amplitude 
exceed  those  of  a  breaking  progressive  wave. 


24 


OlMENSiONUSSS  OISIMNCE 


Fig.  2(a)  —  Blocking  of  a  wave  packet  by  an  advene  noniiniform  cumnt;  the  curvea  diow 
the  aurface  elevation  away  firom  a  mean  level  at  iucceiwve  timaa,  atarting  at  the  bottom. 
Oiatance  ia  raeaaured  in  unita  of  mean  depUi  at  the  left  aide,  ho.  "^e  vertical  diatanre  be* 
tween  profBea  ia  aat  ao  that  a  diaturhance  that  travela  at  a  apeed  >/gho  would  abow  up  ahmg 
a  *  68.1°  angle  in  the  figun.  The  calculationa  wan  run  with  the  unified  wavea  modal,  a  vaii* 
able  breadth  (more  conatricted  toward  the  right)  and  a  ateady  flow  firom  right  to  left  that 
ia  faatar  at  the  right  than  at  the  left.  The  geome^  and  initial  conditiona  are  riiown  in  Fig.  8. 
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Fig.  2(b)  —  Detaib  of  the  aconc  between  x  •  62.6  (grid  point  600) 
end  X  -  137.6  (grid  point  1100)  at  t  -  117.19 
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Fig.  2(e)  —  Detaila  of  the  acene  between  x  -  62.6  (grid  point  600) 
and  X  137.6  (grid  point  1100)  at  t  ■  121.88 
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Fig.  2(e)  —  Detaila  of  the  ecene  between  x  -  62.5  (grid  point  600) 
end  X  -  137.6  (grid  point  1100)  at  t  ••  131.25 
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